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ABSTRACT: We study Hamiltonians consisting of a deterministic term plus a random 
term. Using a diagrammatic approach and introducing the concept of "gluon connected- 
ness," we calculate the density of energy levels for a wide class of probability distributions 
governing the random term, thus generalizing a result obtained recently by Brezin, Hikami, 
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and Zee. The method used here may be applied to a broad class of problems involving 
ON . random matrices. 
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^ , Some four decades ago, Wigner proposed studying the distribution of energy levels of 

a random Hamiltonian given by 

H = cp (1) 
where <p is an iV by N hermitean matrix taken from the distribution 

P{x>) = \e- NtrV M. (2) 
Z 

with Z fixed by J dipP(<p) = 1. This problem has been studied intensively by Dyson, 
Mehta, and others over the years. 2 '^' 4 Two years ago, Brezin and Zee discovered that, 
remarkably, while the density of eigenvalues depends^ on V, the correlation between the 
density of eigenvalues, when suitably scaled, is independent 6 of V. This universality has 
been clarified and extended by other authors/'^ studied numerically, 1 ^ and furthermore, 
shown to hold even when the distribution (2) is generalized to a much broader class of 
distributions. 11 We expect that the discussion to be given below will hold also for this 
broader class of distributions, but for the sake of simplicity we will not work this through 
here. 



In recent work, ' i3 Brezin and Zee have generalized this Wigner problem to the case 
of a Hamiltonian given by the sum of a deterministic term and a random term 

H = H + ip (3) 

Here Hq is a diagonal matrix with diagonal elements q, i = 1,2, ...AT, and <p a random 
matrix taken from the ensemble (2). For the Gaussian case, namely with V(cp) = |v? 2 , 
Pastur 14 has long ago determined the density of eigenvalues. The work described in Ref. 12 
went beyond Pastur's work in that the correlation function between the density of eigen- 
values in the Gaussian case was also determined. More recently, in a work with Brezin and 
Hikami, 13 we managed to determine the density of eigenvalues for V(ip) = \^ + g^ to 
all orders in g. The correlation function was also computed, but only to first order in g. 

This problem of "determinism plus chance" may be regarded as a generic problem 
in physics, and as such represents a significant generalization of Wigner's problem. For 
example, consider an electron moving in a magnetic field and scattering off impurities. 
We note that these "deterministic plus random" problems are considerably more difficult 
than the purely random problems defined in (1) and (2). A standard approach to solving 
the purely random problem involves diagonalizing the random matrix ip and then use 
orthogonal polynomials to disentangle the resulting expression. Clearly, in (3) we cannot 
diagonalize Lp without un-diagonalizing Hq and thus the orthogonal polynomial approach 
fails. 

In this paper, we point out that the problem given in (3) is a special case of a broader 
class of problems involving the addition of random matrices. The deterministic Hamilto- 
nian Hq may in turn be replaced by a random Hamiltonian. Indeed, a deterministic matrix 
is but a special case of a random matrix. We will extend the work of Brezin, Hikami, and 
Zee 13 and determine the density of eigenvalues for the Hamiltonian given in (3) for an 
arbitrary V. 

Our work is inspired by recent advances in the mathematical literature involving the 
theory of non-commutative probability and operator algebra. 15 ' 16 ' 17 A number of physicists 
have already brought these advances to the attention of the physics community. 18 ' 19 ' 20 
While our work is thus inspired, we will not be using the mathematical approach given 
in Ref. 15,16,17, but instead will be based on the diagrammatic approach developed in 
Ref. 12 and subsequent work. 21 ' 13 

Consider a Hamiltonian given by 

H = i Pl + cp 2 (4) 
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with the matrices 2 taken from the probability distribution 

P{<P1,<P2) = ±e- Ntr ^^ +V ^ = Pi(^i)P 2 (^). (5) 

Notice that the probability distribution factorizes. (This is known as "free" in the mathe- 
matical literature.) The problem defined in (3) represents a special case. Previously, with 
D'Anna and with Brezin we have studied the problem given in (4) but for the more difficult 
case 22 ' 23 in which P((fi, 1P2) contains terms linking ipi and if 2- Indeed, a detailed determi- 
nation of the correlation function over all "distance scales" is a non-trivial problem even 
for Gaussian distributions. 23 The discussion in this paper goes through precisely because 
<fl and <f2 do not couple to each other in P(ipi, ¥2)- 

Let us now mention a few necessary definitions. Define the Green's function 

G( *> s (^b) = / /^ifW2)> z _ ( ^ + y2) (6 ) 

The density of eigenvalues is then given by p(fx) = (^tr5(p> — H)"j = —^lmG(fj, + it). In 
this paper we focus on the density of eigenvalues, leaving the correlation for a future work. 
Note that the factors of N are chosen in our definitions such that the interval over which 
p(p) is non-zero is finite (i.e., of order iV°) in the large N limit. 

We may regard the distribution (5) as defining a (0 + 0)-dimensional field theory. The 
Feynman diagram expansion is then simply obtained by expanding G(z) in inverse powers 
of z and doing the integrals in (6). As explained in Ref. 12, it is useful to borrow the 
terminology of large iV quantum chromodynamics 24 from the particle physics literature, 
and speak of quark and gluon lines. See figure (1) for a graphical representation. (It is 
of course not necessary to use this language, and readers not familiar with this language 
can simply think of the diagrams as representing the different terms one encounters in 
doing the integral in (6).) The quark propagator simply comes from the explicit factor 
of z in (6) and is represented by a single line and given by -. The quadratic terms in 
Vi(ipi) and 1^2(932) determine the gluon propagator, represented by double lines. Here, in 
a minor departure from large iV quantum chromodynamics we have two types of gluons, 
corresponding to ip\ and <p2- The gluon propagators are proportional to 

oc S af3 5j6^ (7) 

The non-Gaussian terms in Vi(<pi) and ^(^2) describe the interaction between the gluons. 

The reason that we can solve this problem is because, while the two types of gluons 
have arbitrarily complicated interactions among themselves, they do not interact with each 
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other. Note that while the gluons both interact with the quark, our problem is such that 
we do not have to include quark loops and thus the quark does not induce interaction 
between the two gluons. This is clear from the definition of our problem. Another way of 
saying this is to note that the Green's function may be represented, by using the replica 
trick, as 

G(z)=lim n ^ J Dip^ DipD(pP((p)ip\ipie~ 5^<*=i ^"{z-tptya (g) 

Note that in this language the represent the quark fields and tp the gluon fields. The 
interaction between gluon and quarks are given by ipaPTpa- (Color indices are suppressed 
here.) The interaction of the gluons with each other is determined by P(f)- Since internal 
quark loops are proportional to the number of replicas n, they vanish in the n — > limit. 

Let us then calculate the Green's function, which as usual can be written as (see figure 

2) 

in terms of the one-particle irreducible self energy £* (z) = 5*-E(z). The self-energy is then 
determined by the set of diagrams in figure (3) with the corresponding equation 

=< ^tr<pi > gc + < jjtr<pl > gc G(z)+ < ^tr<p\ > gc G{zf + .... + (1 <-> 2) 

oo 1 

= E < N tr *l >9C G(z) k - 1 + (1 - 2) 

fc=l (10) 

= ^ < ^ tr l^lG>^- 1 ] + ( 1 - 2 ) 

= h [ h G 9 eii h ) - 1] + {1 " 2) 

In order to write this equation, we have to invoke the factorization of P((pi,<P2), which 
tells us that the two kinds of gluons do not interact, and the large iV limit, which tells us 
that the two kinds of gluon lines cannot cross. 

We are led to introduce in (10) the notion of "gluon connectedness," denoted by "gc" 
henceforth. The necessity for this notion is illustrated by the shaded blob describing the 
interaction of the gluons in figure (3d): it should not include the diagram shown in figure 
(4): this class of diagrams is already included in figure (3b). In other words, a gluon 
connected blob with k external gluon lines is such that it cannot be separated into two 
blobs, with k\ gluon lines and gluon lines respectively, (with k\ + k^ = k of course). In 
the last line we have defined the "gluon connected Green's function" 

G gcl (z) =< ^r— L- > gc (11) 



and similarly G gc 2{z). The operations implied in (11) are clearly allowed since < jjtr- > gc 
is a linear operation. Note also that we have not assumed that V a is an even function of 
its argument. In particular, we include a possible tadpole term indicated by < jjtr^p a > gc 
in (10). 

We should emphasize that the shaded blobs include interactions between gluons to all 
orders. It is very complicated, if not hopeless, to calculate these blobs in terms of V\ and 
Vii but fortunately, as we will show below, we do not have to calculate them explicitly. In 
our previous papers, we regarded the cubic, quartic, and so on, terms in V a as interactions 
and proceeded to calculate the Green's function and correlation function in terms of the 
various coupling constants. We follow a different strategy here, and try to express the 
Green's function G(z) directly in terms of G\(z) and G2(z) where 

G a {z) =< hr-^— > (12) 

1\ Z — if a 

for a = 1, 2 are the Green's functions for two separately and purely random problems. (The 
average in (12) is performed with the distribution P a (ipa) = -^-e~ NtrVa ^ Pc ^ Q f course .) i n 
this way, we attempt to bypass having to deal with V\ and Vi altogether. 

To see how to do this, let us go back to the simpler problem defined by (1) and (2). 
Following the same diagrammatic analysis leading to(10) we find that the Green's function 
G(z) and self energy for this simpler problem are related by 

E(z)= ^ [ ^4 )_i] (i3) 

with, evidently, 

G gc (z) =< j^tr-— > gc (14) 

Combining (14) and (9) we find 

^2^(1) = z (15) 
It is convenient to define a "Blue's function" 

B(z) = \G gc (-) (16) 

z z z 

Thus, we learn that the Blue's function is the functional inverse of the Green's function 

B(G(z)) = z (17) 

From the normalization of the probability distribution P(<p) we obtain trivially the "sum 
rule" G(z) — > \ as z — > oo, thus implying that the Blue's function B(z) — > \ as z — > 0. 



Let us now go back to the more involved problem denned by (4). First, we define for 
a = 1, 2 two Blue's functions B a as the functional inverse of G a respectively. We now see 
that (10), when combined with (9), says simply that 

Z + ± = B 1 (G) + B 2 (G) (18) 
Thus, the law of addition for the Blue's function is given by 

B 1+2 (z) = B 1 (z) + B 2 (z)-- (19) 

z 

This equation tells us how to obtain the Blue's function associated with ip\ + ip 2 from the 
Blue's functions associated with (pi and ip 2 . 

The procedure for determining the Green's function and hence the density of eigenval- 
ues of the problem defined by (1) and (2) is then as follows: given the Green's functions 
G\ and G 2 , determine the corresponding Blue's functions B\ and B 2 , calculate B\+ 2 ac- 
cording to (19), then determine the functional inverse of B\+ 2 to find the desired Green's 
function G(z). 

Let us remark briefly on the connection to the mathematical literature. Voiculescu 15 
has introduced the "i?-transform." It turns out that the R function discussed by math- 
ematicians is simply related to B by B(z) = \ + R(z). In fact, we see that the Dyson- 
Schwinger equation (9) when combined with (17) gives simple B{G{z)) = ^|jy + S(z). 
Thus, the R function of the mathematicians is nothing but the self-energy E of the physi- 
cists expressed in terms of different arguments: 

R(G(z)) = E(z) (20) 

Let us now proceed by building up from a few simple examples. In the most trivial 
case, ip is not random at all, but fixed to be a constant c times the unit matrix. Then 
from the Green's function G{z) = we find the Blue's function B(z) = c + \. For a 
slightly less trivial example, let ip be a diagonal matrix with matrix elements given by 
with i = 1, ...,N. Then the corresponding Blue's function is determined by 

1 =z (21) 



N ^ 



N Y B{z) - 6i 

Next, let P((p) be Gaussian (that is, V(ip) = tr^ip 1 ). Then as is well known (see for 
example Ref. 12), the Green's function is determined by 

z = GW + gjjj (22) 
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Substituting z — > B, we obtain immediately that 25 

B(z) = z + - (23) 

z 

Now we are ready to do our first non-trivial problem. Consider the problem defined 
in (3). Since we know from (21) and (23) the Blue's functions corresponding to the two 
terms in the Hamiltonian, we learn immediately from (19) the Blue's function for H: 

B 1+2 {z) = B 1 (z) + z+--- = B 1 (z) + z (24) 

z z 

with B\ determined by (21) with the substitution B — > B\. The desired Green's function 
G(z) is now determined by solving for the functional inverse 

B 1+2 (G) = z (25) 

or equivalently, upon using (24) 

B 1 (G) = z-G (26) 

Let us now evaluate the two sides of this equation with the function Gi(-). We obtain 
immediately 

G{z) = G 1 {z-G{z)) (27) 

precisely the classic result of Pastur which was obtained diagrammatically in Ref. 12. 

After these simple exercises, we can now immediately go on and solve the general 
version of the problem defined in (3): find the density of energy levels of a Hamiltonian 
given by H = Hq + cp with ip drawn from the general distribution (2). With a slight shift 
in notation, let us call the Blue's function associated with Hq and with tp respectively Bq 
and B 2 . Then the Blue's function associated with H is given by B(z) = Bq(z) + B 2 (z) — \. 
Substituting z — > G[z) (where G[z) is the unknown Green's function associated with H), 
we find immediately that 

B (G) = z + 1 - B 2 (G) (28) 
Anticipating the next step, we define 

E(z) = B 2 (G(z)) - (29) 

Let us now evaluate both sides of (28) with the Green's function Gq(-) associated with Hq. 
We find instantly that 

GW = Go( ,-E W ) = l£— L_ (30) 



7 



Let us repeat this trick: rewrite (29) as -B 2 (G(z)) = E(z) + and evaluate both sides 
with the function G^O- We obtain 

G(z) = G 2 (E(z) + ^y) (31) 

These two equations, (30) and (31), allow us to determine the two unknown functions 
G(z) and provided we know the Green's function G 2 (-2)- But what is the Green's 

function ^2(2)? It is just the Green's function associated with the random matrix if 
drawn from the general distribution (2). But this was obtained by Brezin et al 5 almost 
twenty years ago. These authors told us that (for V{z) an even polynomial for the sake of 
notational simplicity) 

G 2 (z) = \[V'{z) - P(z)Vz 2 -a 2 ] (32) 

Here V'(z) = %, P(z) is a polynomial, and a determines the endpoints of the spectrum of 
eigenvalues. The quantities P(z) and a are determined 26 by the "sum rule" that G2 — > \ 
as z — > 00. 

Thus, in summary, for any distribution defined by V(cp), we can determine the Green's 
function G{z) and hence the density of eigenvalues by solving (30) and (31). It is clearly 
convenient to define a(z) = + t^t-t- We can then simplify (31) slightly to 27 

P 2 (a) (a 2 - a 2 ) = (V'(a) - 2G) 2 (33) 

Thus, we can use (33) to determine cr, and hence S, in terms of G. Plugging this into (30) 
then gives us an equation for G. 

As mentioned earlier, Brezin, Hikami, and Zee 13 recently used the equation of motion 
method and a detailed diagrammatic analysis to determine the Green's function for the 
problem in (3) with the distribution defined by V(<f) = ^ip 2 + gip^. It is straightforward, 
although slightly tedious, to verify that for this simple case, (33) reduces to equation (4.15) 
in Ref. 13. The analysis given here is considerably simpler. 

In practice, with an arbitrary set of ej, (30) can only be solved numerically. One 
relatively simple example involves taking half of the q's to be equal to +e and the other 
half to be equal to — e. Thus, before the introduction of the random term ip into the 
Hamiltonian, the spectrum consists of two levels, which may for example represent the 
two lowest Landau levels, corresponding to spin up and spin down, in a spin-dependent 
quantum Hall fluid. 28 ' 29 ' 30 . The randomness will then broaden the two levels. 
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Clearly, by repeating the discussion given here, we can add an arbitrary number of 
random Hamiltonians together. The procedure is defined by (19). The deterministic plus 
random Hamiltonian studied here is just a special case. It will be interesting to see if 
this work can be generalized to a study of the universal correlation function discussed in 
Ref. 6,12,13. 
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FIGURE CAPTIONS 

Fig 1. Feynman rules: (a) quark propagator, (b) gluon propagator, (c) quark gluon 
vertex, (d) gluon interaction, illustrated here with a gtp^ vertex. 

Fig 2. Quark propagator and one-particle-irreducible self energy. 

Fig 3. Quark self energy: the gluons shown explicitly are all of type 1. There are of 
course also type 2 gluons inside the quark propagator G. 

Fig 4. A class of diagrams not included in (3d). 
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